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In the standard model (SM), the electromagnetic ra- 
diative decays of the b quark, b — >• s'y and b — > dj, pro- 
ceed via a loop diagram at leading order. A wide variety 
of new physics (NP) scenarios such as supersymmetry 
may cause new contributions to the loop [l|48| at the 
same order as the SM, resulting in significant deviations 
for both the branching fractions and the direct CP asym- 
metry 



A C p = 



T[b -> (s + d) 7 ] - T[b -> (a + d)j] 



T[b -> (s + d) 7 ] + T[b -> (S + d) 7 ] 

Inclusive hadronic branching fractions (BF) B(B 
X s ^) and B(B X^j) can be equated with the pertur- 
batively calculable partonic BF B(b — > sj) and B(b — > 
dj) at the level of a few percent [9|, allowing theoreti- 
cally clean predictions. At next-to-next-to-leading-order 
(four-loop), the SM calculation gives B(B — > X s j) — 



(3.15 ±0.23) x 1(T 4 (£ 7 > 1.6 GeV) 



where E~ is 



B meson. B(B —> Xdj) is suppressed by a factor of 
|Vt<i/Vt s | 2 ~ 0.04, where Vij are the elements of the 
Cabbibo-Kobayashi-Mashawa(CKM) quark-mixing ma- 
trix. NP with non-minimal flavor violation can also sig- 
nificantly enhance Aqp [lit , which is approximately 10~ 6 
in the SM [l2l - fl4l |. Consequently the precision measure- 
ment of these decays has long been identified as impor- 
tant in the search for NP. They are central to the program 
of the future Super B factories ll|-I7|, which will probe 
NP mass scales up to 100 TeV. 

In this letter, new precise measurements of B(B — > 
X s -f) and Aqp are presented. The analysis has been sig- 
nificantly improved from our previous result [l8j |. which 
it supersedes. In addition, the shape of the photon en- 
ergy spectrum is measured in the B meson rest frame. It 
is insensitive to NP [l!| but can be used to determine the 
Heavy Quark Expansion parameters mt, and 2^, 21 



the photon energy measured in the rest frame of the 



related to the mass and momentum of the b quark within 
the B meson. These parameters arc used to reduce the 
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uncertainty in the extraction of the CKM elements \V c b\ 
and | V u b from semi-leptonic B meson decays [22h25|. 

This Letter summarizes a fully inclusive analysis of 
B ->• X s j decays collected from e + e~ T(45) -> BB 
events. Full details are given in Ref. (2(|. The photon 
from the decay of one B meson is measured, but X s is 
not reconstructed. This avoids large uncertainties from 
the modeling of the X s system, at the cost of large back- 
grounds, which need to be strongly suppressed. The prin- 
cipal backgrounds are from other BB decays containing 
a high-energy photon and from continuum qq (q = udsc) 



and 



events. The continuum background, including 



a contribution from initial-state radiation, is suppressed 
principally by requiring a high-momentum charged lep- 
ton ("lepton tag") from the non-signal B decay, and 
also by discriminating against events with a more jet-like 
topology. The BB background to high-energy photons, 
dominated by ir° and rj decays, is reduced by vetoing 
reconstructed 7r° or r\ mesons. The residual continuum 
background is subtracted using off-resonance data col- 
lecdtcd at a center-of-mass (CM) energy 40MeV below 
the T(4S f ), while the remaining BB background is es- 
timated using a Monte Carlo (MC) simulation that has 
been corrected using data control samples. The photon 
energy spectrum is measured in the T(AS) rest frame. 
Quantities measured in this frame are denoted by an as- 
terisk, e.g., E*. 

The data were collected with the BABAR detector [53] 
at the PEP-II asymmetric-energy e + e~ collider. The 
on-resonance integrated luminosity is 347.1 fb _1 , corre- 
sponding to (382.8 ± 4.2) x 10 6 BB events. Addition- 
ally, 36.4 fb -1 of off-resonance data are used. The BABAR 
MC simulation, based on GEANT4 [Hj|, EVTGEN 
and JETSET 30], is used to generate samples of B + B~ 
and B°B° (excluding signal channels), qq , t + t~, and 
signal events. The signal models used to compute effi- 
ciencies are based on QCD calculations in the "kinetic 
scheme" [13], "shape function scheme" 21], and in 
an earlier model These calculations approximate 

the X s resonance structure with a smooth distribution 
in the hadronic mass mx a ■ The portion of the mx a 
spectrum below l.lGeV/c 2 , where the if* (892) domi- 
nates, is replaced by a Breit-Wigner if*(892) distribu- 
tion. The analysis is performed "blind" in the range 
1.8 < E* < 2.9 GcV; that is, the on-resonance data are 
not examined until all selection requirements are finalized 
and the corrected BB backgrounds determined. The sig- 
nal range is limited by large BB backgrounds at low E*. 

The event selection begins by requiring at least one 
photon candidate with 1.53 < E* < 3.50 GeV. A photon 
candidate is an electromagnetic calorimeter (EMC) en- 
ergy cluster with a lateral profile consistent with that of 
a single photon, isolated by 25 cm from any other clus- 
ter, and well contained in the calorimeter. Photons that 
are consistent with originating from an identifiable ir° or 
rj — > 77 decay are vetoed. Hadronic events are selected 



by requiring at least three reconstructed charged parti- 
cles and the normalized second Fox- Wolfram moment Ri, 
to be less than 0.9. To reduce radiative Bhabha and 
two-photon backgrounds, the number of charged parti- 
cles plus half the number of photons with energy above 
0.08 GeV is required to be at least 4.5. 

About 20% of B mesons decay semileptonically to ei- 
ther e or fi. Leptons from these decays are emitted 
isotropically and tend to have higher momentum than 
the continuum background in which the lepton and pho- 
ton candidates also tend to be anti-collinear. A tag- 
ging lepton (£ — e, /i) is required to have momentum 
p* > 1.05GeV/c and an angle relative to the photon 
cos^ > —0.7 to suppress this background. It does 
not compromise the inclusiveness of the B — > X s j se- 
lection since it comes from the recoiling B meson. The 
presence of a relatively high-energy neutrino in semilep- 
tonic B decays is used to further suppress the background 
by requiring the missing energy of the event to satisfy 
£* iss >0.7GeV. 

The sample is separated into electron and muon tags. 
For each, p* and cos9* e are then combined in a neu- 
ral network (NN) with eight event-shape variables that 
exploit the difference in topology between isotropic BB 
events and jet-like continuum events. The NN is trained 
to separate signal-like events from continuum background 
using MC samples. The BB background sample is ex- 
cluded from the training because it is used for back- 
ground subtraction and is topologically similar to the 
signal. 

The selection criteria are optimized for statistical pre- 
cision. This was done iteratively for five variables: the 
two NN outputs, the energies of the lower-energy pho- 
ton in the n° and r\ vetoes, and i^iss- The signal effi- 
ciency for the entire selection depends on E*, falling at 
lower values. This effect is significantly reduced from our 
previous analysis, lessening the uncertainty due to the 
assumed signal model ("model-dependence"). The effi- 
ciency integrated over the range 1.8 < E* < 2.8 GeV is 
about 2.5%, while only 0.0005% of the continuum and 
0.013% of the BB background remain in the sample. 

The remaining continuum background is estimated 
with off-resonance data scaled to the on-resonance lu- 
minosity and adjusted to account for the 40 MeV CM en- 
ergy difference. The BB background is estimated with 
the BB MC sample. It consists predominantly of pho- 
tons originating from 7r° or rj decays (~ 80% in the signal 
region), electrons (« 10%) that are misreconstructed, not 
identified, or undergo hard bremsstrahhmg, uj and r/ de- 
cays (« 4%), and n's (« 2%) that fake photons by anni- 
hilating in the EMC. Each of the significant components 
is corrected by comparison with data control samples. 

The 7T° and rj background simulations are compared to 
data using the same selection criteria as for B — > X s j 
but removing the ir° and rj vetos. For this compar- 
ison the high-energy photon requirement is relaxed to 
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E* > 1.03 GeV to increase the size of the sample. The 
yields of tt° and r\ are measured in bins of E* ^ by 
fitting the 77 mass distributions in on-resonance data, 
off-resonance data, and BB simulation. Correction fac- 
tors to the 7T° and r\ components of the BB simulation 
are derived from these yields. An additional correction 
is applied to account for data-MC differences in the low- 
energy photon detection efficiency This has an opposite 
effect on the control-sample 7r° and 77 selection than on 
the standard event selection, where finding a it or 77 re- 
sults in the event being vetoed. 

As an antineutron control sample could not be isolated, 
this source of BB background is corrected by comparing 
in data and simulation the inclusive antiproton yields in 
B decay and the EMC response to p's, using A — > pn + 
samples. The misreconstructed electron background is 
measured using B — > X J/ip(e + e~) data. This sample 
closely models the particle multiplicity in B — > X s "f 
events. Bremsstrahlung in the detector is reliably simu- 
lated by GEANT4, so no correction is necessary. The 
small contributions from u> and rj decays are corrected 
in bins of E* using inclusive B decay data. Nearly all of 
the tagging leptons arise from B — > X c lv. The yield of 
such events in the simulation is corrected as a function of 
lepton momentum according to previous BAB4R measure- 
ments 31,33]. The complete BB background estimation 
incorporates the correction factors and uncertainties and 
includes correlations between E* bins. The dominant un- 
certainties originate from the 7r°, 77, and misreconstructed 
electron corrections. 
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FIG. 1: The measured Ej photon energy spectrum after back- 
ground subtraction, uncorrected for efficiency and resolution 
smearing. The inner error bars are statistical only, while the 
outer include systematic errors added in quadrature. 

Figure Q] shows the measured E* spectrum after sub- 
tracting both continuum and BB backgrounds. The sys- 
tematic errors are due to the BB subtraction uncertainty. 
The region 1.53 < E* < 1.80 GeV is dominated by BB 
background, while the higher-energy range 2.9 < E* < 



3.5 GeV contains only continuum background. These 
regions are used to validate the background subtrac- 
tion procedure. In the higher-energy range there are 
— 100± 138(stat) events. In the lower-energy region there 
are 1252 ± 272(stat) ± 841(syst) events. Allowing for an 
average of 275 signal events from a range of plausible 
signal models, and for correlations between the bins, the 
latter result is consistent with zero to within one stan- 
dard deviation (ler). 

To extract BF's and the shape of the spectrum, it is 
necessary to first correct for efficiency. Theoretical pre- 
dictions are made for the true B 7 in the B meson rest 
frame, whereas the E* is measured in the T(AS) frame. 
Hence it is also necessary to correct for the asymmet- 
ric EMC resolution and the Doppler smearing due to 
the motion of the B meson in the T(4S) rest frame. 
The efficiency and smearing corrections depend upon the 
assumed signal shape. In both the kinetic and shape 
function schemes this shape is parameterized by mb and 
/4- The Heavy Flavor Averaging Group (HFAG) [H] 
has extracted values and uncertainties in the kinetic 
scheme by fitting moments of inclusive distributions in 
B — > X c lv decays and previous B —y X s ~/ measure- 
ments, and has also translated them to the shape function 
scheme. These results define the nominal signal model 
(kinetic scheme) used for the BF measurement, along 
with a mo del- dependence uncertainty (kinetic and shape 
function schemes). To provide an independent measure- 
ment of the shape of the spectrum, the measured spec- 
trum is unfolded using an iterative technique that reduces 
sensitivity to the signal model. In this case the initial sig- 
nal model and model-dependence uncertainty are based 
on the data rather than the HFAG parameters. The ef- 
fects of efficiency and smearing cancel in the Aqp mea- 
surement so it is extracted directly from the measured 
E* yield separated by lepton tag charge. 

The BF is computed from 



B(B -> X s+dl ) 



aS/(2N B ^e s 



5 ), 



where S is the signal yield integrated over the E* ranges 
1.8,1.9,2.0 to 2.8 GeV, e s i g is the signal efficiency, and 
N bb i s the number of BB pairs in the sample. The 
factor a, which is close to unity, corrects for resolution 
and Doppler smearing and is computed with the nominal 
signal model. The model-dependence errors on the BF 
associated with the efficiency and the smearing correc- 
tion are fully correlated. The results for the three energy 
ranges are given in Table „ The BF's have been corrected 
by a factor 1/(1 + {\V t d\l\V ts |) 2 ) = 0.958 ± 0.003 to 
remove the contribution from b — > dj. The most signifi- 
cant systematic error is from the corrections to the BB 
background simulation, which in the range 1.8 GeV < 
E~f < 2.8 GeV contributes 7.8% to a total systematic 
uncertainty of 9.0%. Additional contributions added in 
quadrature, all energy-independent, arise from uncer- 
tainties in the selection efficiency (3.1%), predominantly 
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TABLE I: The measured BF, first, and second moments (±stat ± syst ± model) for different ranges of E 1 in the B rest frame. 
Correlations between the energy ranges are given in Ref. [2(|. 



E 1 Range ( GeV) 


B(B -»• X s j) (10" 4 ) 


(E- f ) (GeV) 


((E, - (E,) f) (GeV 2 ) 


f .8 to 2.8 


3.21 ±0.f5± 0.29 ±0.08 


2.267 ± 0.019 ± 0.032 ± 0.003 


0.0484 ± 0.0053 ± 0.0077 ± 0.0005 


f .9 to 2.8 


3.00 ±0.14 ±0.19 ±0.06 


2.304 ± 0.014 ± 0.017 ± 0.004 


0.0362 ± 0.0033 ± 0.0033 ± 0.0005 


2.0 to 2.8 


2.80 ±0.12 ±0.14 ±0.04 


2.342 ± 0.010 ± 0.008 ± 0.005 


0.0251 ± 0.0021 ± 0.0013 ± 0.0009 



due to the high-energy photon and NN selections, the 
scmileptonic BF for B meson decays, and the modeling of 
the X s system. Correlations between the BB and the sig- 
nal efficiency systematic errors contribute an additional 
2.9% uncertainty. Finally, there is a 1.1% uncertainty in 



N 



BB ■ 




Ey (GeV) 

FIG. 2: The E^ photon energy spectrum corrected for effi- 
ciency, resolution, and Doppler smearing, shown as a partial 
branching fraction AB. The inner error bars are statistical 
and the outer include systematic errors added in quadrature. 
The vertical line shows the boundary between the lower con- 
trol region and the signal region. The curve is the kinetic 
scheme model using HFAG world average parameters, nor- 
malized to data in the range 1.8 < E? f < 2.8 GeV. 

To obtain an E 1 spectrum in the B rest frame, the E* 
spectrum shown in Fig. Q] is corrected for selection effi- 
ciency, and the resolution smearing and Doppler smear- 
ing are unfolded. A simplified version [3o| of an itera- 
tive unfolding technique [36j is used. The method starts 
with an initial signal model that, when passed through 
the detector simulation and event selection, closely re- 
sembles the data (shape function scheme with nit, — 
4.51 GeV, [j% = 0.46 GeV 2 ). This model is used to cor- 
rect for efficiency and unfold the data. A fraction, de- 
termined by a bin-dependent regularization function, of 
the difference between the unfolded data and the initial 
signal model is used to adjust the signal model, and the 
process is iterated until it converges. Only one itera- 
tion is necessary. The results are shown in Fig. [5J This 
technique preserves fluctuations in the spectrum and re- 



duces the model error. The model dependence uncer- 
tainty is computed using an initial model that is approx- 
imately lcr lower than the data in Fig. [T] in the region 
with significant BB background (1.8 < E* < 2.1 GeV). 
The error is the absolute value of the difference bin by 
bin after unfolding. It is small except near the kine- 
matic limit, _E 7 w ms/2, where the sharply falling edge 
leads to strongly anti-correlated differences in adjacent 
bins. To reduce this effect, the 100-MeV bins between 
2.4 and 2.8 GeV are combined into 200-MeV bins. The 
spectral shape and the full covariance matrix, provided 
in Ref. |26j . are used to compute the first and second 
moments in Table „ They can also be used to fit any 
theoretical prediction for the spectral shape. The BF's 
computed from the sum of the AB in Fig.[2]are consistent 
with the values given in Table „ [26| . 

Finally the E* sample is divided into B and B decays, 
using the charge of the lepton tag, to measure A™p as (B — > 
X s+dl ) = (N+ - N-)/(N+ ± N~) where AT+(") are 
the positively (negatively) tagged signal yields. Aqp is 
then given by Aqp = v4™p as /(l — 2ui) where ui is the 
mistag fraction. To maximize the statistical precision 
a requirement of 2.1 < E* < 2.8 GeV is made. This 
is determined from simulation and does not bias the 
SM prediction for the asymmetry 
N+ = 2620 ± 158(stat) and N" 



37J. The yields are 
= 2389 ± 151(stat). 
The bias on Aqp due to charge asymmetry in the de- 
tector response or BB background is measured to be 
AAg£ as (£_-> X s+dl ) = -0.004 ± 0.013, using events 
in the BB control region to check for a background 
asymmetry, and using several event samples (e + e~ — > 
e+e~7, e+e" -)• hhj and B ->■ J /ip (£+£-)) to 
check for a lepton tag asymmetry. The mistag frac- 
tion ui — 0.133 ± 0.006 is dominated by B°B° mixing, 
which contributes 0.093 ± 0.001 34], with an additional 
0.040 ± 0.005 arising from wrong-sign leptons from the 
B decay chain and from misidentifcation of hadrons as 
leptons. After correcting for charge bias and mistagging 
it is found 

A CP = 0.057 ± 0.060(stat) ± 0.018(syst). 

The systematic error includes relative uncertainties from 
the BB background subtraction (2.2%) and mistagging 
(1.8%). The uncertainty due to differences in the B — > 
X s -f and B — > Xjy spectra is negligible. 
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In summary, the photon spectrum of B — ¥ Xg+d'y de- 
cays has been measured and used to extract the branch- 
ing fraction, spectral moments, and Aqp- Previous in- 
clusive measurements of B — > Xsj have been presented 
by the CLEO [H, BaBAR [3, and Belle Col- 
laborations. The measured branching fraction B(B — > 
X sl ) = (3.21 ± 0.15 ± 0.29 ± 0.08) x 10~ 4 (1.8 < E 1 < 
2.8 GeV) is comparable in precision to the Belle result, 
(3.36 ± 0.13 ± 0.25 ± 0.01) x 10~ 4 , but with a dataset 
that has 60% smaller integrated luminosity. The BF 
for 1.8 < E 1 < 2.8 GeV is extrapolated to the range 
£ 7 > 1.6 GeV using a factor of 1/(0.968 ± 0.006) de- 
termined by HFAG. This results in B(B X s j) = 
(3.31 ± 0.16 ± 0.30 ± 0.09) x 10~ 4 for £ 7 > 1.6 GeV, 
in good agreement with the SM prediction. The extrap- 
olated B(B — > X s j) can be used to constrain NP. For ex- 
ample, in a type-II two-Higgs-doublet model HI the 
region M H ± < 327 GeV is excluded independent of tan/3 
at 95% confidence level. This limit is far more stringent 
than that from direct searches at the LHC 4l|, |42| . The 



Acp measurement is the most precise to date and can 
be used to constrain non-minimal flavor-violating mod- 
els 11|. The measured moments and spectra provide 
input to improve the precision on the HFAG estimation 
of mi, and /xj, which will result in a reduced error on 
\V u b\- Finally, the improved technique presented in this 
paper can be applied with increased precision at future 
Super B factories. 
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